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We introduce a hadronic scale i?o through the force F{r) between static quarks at intermediate 
distances r. The definition F(Ro)Rq = 1.65 amounts to Rq ~ 0.5 fm in phenomenological potential 
models. Since Rq is well defined and can be calculated accurately in a Monte Carlo simulation, it is 
an ideal quantity to set the scale. In SU(2) pure gauge theory, we use new data (and Rq to set the 
scale) to extrapolate F(r) to the continuum limit for distances r = 0.18 fm to r = 1.1 fm. 
Through Rq we determine the energy scale in the recently calculated running coupling, which used 
the recursive finite size technique to reach large energy scales. Also in this case, the lattice data can 
be extrapolated to the continuum limit. The use of one loop Symanzik improvement is seen to reduce 
the lattice spacing dependence significantly. 

1. Introduction 

In the hmit when the mass of a quark h becomes large compared to typical QCD- 
scales (such as the proton's mass), bound states hh are expected to be described by an 
effective nonrelativistic Schrodinger equation[l]. The nonrelativistic potential is given 
by the energy of static fundamental charges separated a distance r. In the real world, 
such a description should be approximately valid for the bb and maybe even for the cc 
spectra. 

In fact, both these spectra can be described by one effective potential. Since (as one 
knows a posteriori from a successful model) the bb states have rms-radii of 0.2 fm to 
0.7 fm and the cc states have a size of 0.4 fm to 1.0 fm, the spectra determine the 
effective potential in the range of r ~ 0.2 fm to r ~ 1 fm. Although we must remember 
that the relationship between the static QCD potential and the effective potential used 
in phenomenology is not well understood, this discussion suggests that the distance 
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range where we have the best information on the force F[r) between static quarks is 
at distances of around 0.5 fm. 

In lattice gauge theory calculations, we need to fix one dimensionful quantity in order 
to set the overall scale. In the pure gauge theory, this is usually done through the 
string tension, defined as K = linir^oo F[r). The limiting procedure is not easy to do 
because the statistical and systematic errors on the force increase with the distance. It 
is therefore superior to use the force at intermediate distances to determine the scale: 
one calculates R{c) such that F[R[c))R[cy = c. Because of the above mentioned 
arguments, we chose c = 1.65, which corresponds to Rq = i?(1.65) = 0.49 fm in the 
Cornell [2] and the Richardson [3] model and very similar values in the other successfuU 
models. 

The statistical and the systematic errors of Rq are quite moderate. This opens the 
possibility of extrapolations to the continuum limit. 

In section 3, we determine the force in the continuum limit of the SU(2) Yang Mills 
theory: the dimensionless quantity H[x) = F[r)r'^\r=xRo is calculated at five different 
values of the lattice spacing and then extrapolated to lattice spacing zero. It is well 
determined just about in the same range as the phenomenological potential. We 
compare H[x) to the phenomenological potential, the short distance perturbative 
form and the expected large distance form. The latter is given by an effective bosonic 
string theory including its universal leading corrections. 

Our main motivation, however, was to set the energy scale in the calculation of the 
running coupling. In ref. [5] a physical coupling was defined for the pure gauge theory 
in finite volume, that depends only on one scale which is the size of the box L. One 
can calculate the step scaling function which gives the change of the coupling g^{L) 
when the scale L is changed by a factor s. This computation was done in ref. [6] and 
the lattice results could be extrapolated to the continuum. As a result g^{L) is known 
for LjU = 1, 0.500(23), 0.249(19), 0.124(13) and 0.070(8). Here U is defined by 
g'^{Ls) = 4.765. 

The last step is to calculate the relation of Lg, to some physical scale. To this end we 
previously[6] used the string tension but we could only give a crude estimate of Lg^/K. 
In section 4 we will instead calculate Lg/ Rq and extrapolate it to the continuum. In the 
part of the calculation that involves g^[L)^ we used the Symanzik improved action[5] 
where the 0{a/L) lattice artifacts are removed to one loop order. Comparing this to 
the results without improvement, we observe a very significant reduction of the lattice 
artifacts. 
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2. Rq — a. Precise Low Energy Scale 



In order to set the scale in a lattice gauge theory calculation, we need to determine a 
physical dimensionful quantity. As to today's knowledge, the force F[r) between static 
quarks is the quantity which can be calculated most precisely. In order to avoid an 
extrapolation in the distance r, we define a hadronic length-scale R{c) which depends 
on a dimensionless parameter c through 

F{r)l=nic) = c. (2.1) 

As discussed in the introduction, choosing e.g. c = 1.65 corresponds to 

i?(1.65) = i?o ^ 0.5 fm . (2.2) 

The advantages of this choice are: 

• R{c) is defined precisely, both in the pure gauge theory and in the theory with 
dynamical fermions. In particular, it does not refer to the quenched approxima- 
tion. 

• In contrast to the string tension it does not contain residual systematic errors 
that originate from assumptions on subleading terms in the force at intermediate 
distances. 

• It can be calculated with good statistical precision. Therefore its use to set 
the scale for other physical observables will not introduce significant additional 
uncertainties into the observables of interest. This is of special importance when 
one wants to extrapolate lattice results to the continuum. Examples are given 
in sections 3.2 and 4. 



2.1 Relationship to other Observables 

As in the case of the string tension, one has to use the effective potential in order to 
relate Rq (or more generally R{c)) to experimental observables. Using the Richardson 
potential, we obtain Rq = 0.49 fm. So e.g. the relationship to the IP - IS splitting in 
the cc system Am and to the proton's mass nip is given by 

Am Ro = 1.14, Am i?(l) = 0.72 (2.3) 
nip Ro = 2.33, nip R{1) = 1.47 (2.4) 

Very similar numbers are obtained when one uses other QCD-inspired phenomenolog- 
ical potentials, whereas the logarithmic potential gives values that are lower by about 
10%. 
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Once one can calculate in full QCD it will be interesting to check eq.'s (2.3,2.4) and 
one can get a certain measure of the precision of the nonrelativistic description. 



2.2 Numerical Determination of R{c) 



It is of vital importance that the quantity that one uses as a reference scale can be 
computed accurately. We demonstrate that this is the case for R{c) in the following. 
Three steps have to be discussed in the computation of R{c) in lattice units: 1) the 
calculation of the potential at a certain distance r/a, 2) the definition of the force 
through a finite difference from the potential and 3) the interpolation of the force 
which is necessary to determine R[c)/a through eq. (2.1). 

Step I) is important but rather technical. It has already been discussed in detail in 
[7]. Our specific method is described in the appendix. Here, we summarize only the 
conclusion: using a variational principle and a large enough range of t for the smeared 
Wilson loops, one obtains an upper bound Kt(r) and a lower bound Vi{f) for the 
potential. We point out that it is advantageous to have larger values of t than the 
ones used in [7], since then Kt(r) and Vi{f) are close to each other and have moderate 
statistical errors. 

Step 2): given the potential, the force at distance r/ and along the orientation d is 
computed from 



Fdin) = \d\-'mf)-Vi{f-d)], (2.5) 

n = [ iir \d\-^ {G{f) - G{f - d)) ]-^^' (2.6) 

Gir) - a-' r n,tiCos(r,%/a) 

^ ^ " 1. (27r)3 4E,iisin2(%/2) ^'"^^ 

Here, the argument r/ is chosen such that Fj^rj) is a treelevel improved observable, 
i.e. to order (/q we have Fj(^ri) = (/Q/(47rr|). This definition removes the (9((a/r)^) 
lattice artifacts that would be present in the naive choice r = \ f — d/2\ as argument 
of F. No fit is necessary to achieve this. Note in addition, that the exact choice of 
the argument is irrelevant when the force becomes constant. Eq.(2.5) is a convenient 
definition that greatly reduces the lattice artifacts in the force. We stress, however, 
that the remaining lattice artifacts have to be controlled by extrapolating simulation 
results with various (small) a to the continuum. This will be described in the following 
sections. 

We calculated the central values of the force from Vi. In order to cover the systematic 
error of the force that is due to using eq. (A7) at finite values of t, we always repeated 
the calculation of the force with Vi{f) — )• Kt(r) and added the resulting difference of the 
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central values to the largest statistical errors of the force. The latter originated from 
Vi{f). Although there is no reason that inserting Kt(r) instead of Vi{f) into eq. 2.5 will 
always result in a larger value for the force, this is generally the case in the data. This 
reflects mainly the fact that the gap eq. (A5) decreases with growing \ f\ and therefore 
Kt(^) — V;(r) increases with |r|. 

Step 3) consists of interpolating the force. Here one can simply use the two neighboring 
points and the interpolation function F[r) = /i + /2 . As this is locally an 
excellent approximation to the r-dependence, the interpolation error is small: it can 
be estimated from the change that arises from adding a term /s r""* and taking a third 
point. In the applications discussed below, this change was added as a systematic 
error. It is well below the statistical uncertainty. 



3. Reconstruction of the Force in the Continuum 

As a first application of setting the scale through i?o, we describe here the calculation 
of the force in the continuum limit. 



3.1 The Monte Carlo Simulation 

In order to be able to extrapolate to the continuum limit, we performed computa- 
tions at different values of the lattice spacing. We simulated L"^ lattices with periodic 
boundary conditions and the standard Wilson action [8]. According to the state of the 
art, we used a hybrid algorithm performing N exactly microcanonical overrelaxation 
sweeps followed by one Creutz heatbath sweep. Table 1 summarizes the parameters 



/3 


L/a 




tmax / ^ 


M 


N 


F 


sweeps 


2.50 


16 


5, 15, 40 


8 


4 


4 


1 


85k 


2.55 


20 


4, 30 


9 


6 


5 


2 


30k 


2.60 


24 


5, 40 


10 


7 


6 


3 


56k 


2.70 


32 


3, 16, 40 


8 


8 


8 


2 


27k 



Table 1: Parameters of our simulations. Measurements were taken every F x (N + 1) sweeps. In 
order to reduce the auto correlations between measurements, we cycled through the eight sublattices 
A J (cf. the appendix) from measurement to measurement. Si gives the number of smearing iterations 
as discussed in the appendix, tmax denotes the maximum time-extent of the Wilson loops and the 
maximum spatial separation is rmax = M2\/3a. 

used in the simulations. The largest system was simulated on the CERN-IBM, where 
a speed of about 65 Mfiop/s could be achieved. 
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Integrated autocorrelation times were estimated for all smeared loops using jacknife 
binning. They are e.g. roughly 10 sweeps at /3 = 2.5 and 14 sweeps at /3 = 2.6 for the 
smallest smeared Wilson loops and decrease with increasing loop size. 

In table 3, we give the results for the force on all our four lattices. The errors quoted 
are the jacknife errors[9] on eq. (2.5) plus the change in the central value of the force 
when we replace Vi{f) — )• Kt(r). For /3 = 2.85, the corresponding numbers were read 
off from table 4 of ref. [7]. A number of our /3 = 2.70 results can be compared with the 
ones in [10]. Our central values are all smaller. This difference is somewhat outside 
of the error bars. Although it is possible that the largest difference is only a typing 
mistake in [10], it appears that the general difference is due to the smaller range in t 
used in [10]. The relatively large error bars at /3 = 2.85 reflect the fact that the gap in 
lattice units is small. Hence Kt(r) — Vi{f) becomes significant at such a small lattice 
spacing and the calculation of the force is difficult at /3 = 2.85. 

3.2 Taking the Continuum Limit 

Once we have set the physical scale through i?o, the force is entirely described by a 
dimensionless function: 

H{x) = F{r)r^\r=,R,. (3.1) 

At different values of the lattice spacing and different orientations di , we obtain the 
lattice approximations Hj;[x^a) to H[x). From Symanzik's discussion of the cutoff 
dependence of loop integrals [11], one expects that the lattice approximations converge 
to the continuum with corrections that are roughly proportional to a^: 

H^-^{x,a/Ro) = H{x) + 0{a''/Rl) (3.2) 

The lattice results are plotted in fig. 1 for five sample values of x as a function of the 
square lattice spacing [a/i?o]^- The extrapolations according to eq. (3.2) are shown as 
well. Within our precision, no significant a-dependence is seen and the extrapolations 
to the continuum are stable and clearly do not depend on the functional form of 
the assumed a-dependence. To a large extent this is due to the tree level improved 
definition of the force (cf. section 2.2). 

Fxtrapolations were done for a range of 0.35 < x < 2.1. At smaller values of x, we 
have not performed the extrapolation since there we have only two points in a. For 
X > 1.1 there is data with /3 < 2.7 only. 
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Figure 1: Continuum extrapolation for H(x) with x = 1.6 (A) x = 1.4 (□), and x = 
0.8, 0.6, ,0.4, 0.35 (o) from top to bottom. The circles are for the force along orientation di, 
the squares correspond to orientation d2 and the triangles to orientation d^. 

3.3 Discussion of the Force 

The continuum results for H[x) are plotted in fig. 2. It is determined accurately in 
about the same distance range, where the universality of the phenomenological force 
holds [4]. 

We compare this SU(2) Yang-Mills force with the phenomenological model by Richard- 
son[3]. As can be seen e.g. in fig. 1 ref. [4], other phenomenological models do not 
differ much from the Richardson model in the range of x discussed here. We note, 
however, that the force of the Martin[I2] and the logarithmic potential [13] differs 
from the Richardson force as much as the SU(2) Yang-Mills force. In this sense, even 
without fermions and for gauge group SU(2), the static force compares quite well with 
the phenomenological forces. 
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Figure 2: H(x) = _F(r)r^ |r=j;_Ro as a function of x. Symbols as in fig. 1. The Richardson model is 
shown as full line and the bosonic string model (normalized by -ff(l) = 1.65) as dashed curve. 

In the figure, we also show the prediction of the bosonic string model normalized at 
X = 1. It includes besides the string tension the universal 7r/(I2r^) correction to the 
force. This correction should be there at large r [14]. As already noted earlier [15], this 
form gives an excellent effective representation of the force also for rather small values 
of X. Note that previous claims in the literature[16], that the Yang-Mills Wilson loops 
are effectively described by a fermionic string model are based on the subtraction of 
perturbative contributions. 

3.4 The Running Coupling agg(r) 

The quantity H[x) may be used to define a physical running coupling 

4 

'=^qq{^') = '^H{x), r = xRq. (3.3) 
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It is of interest to see, in how far we have reached the perturbative region in x. To 
check this, we start with our value of agg(r) at the smallest value of r, and integrate 
the 1-loop and the 2-loop perturbative renormalisation group equations towards larger 
r. This is shown in fig. 3. We see, that the 2-loop /3-function describes the evolution 
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Figure 3: agg(f) together with 1-loop (dashed) and 2-loop (dotted-dashed) renormalisation group 
evolution starting at the smallest value of the coupling. 

of the coupling for only a change in scale up to a factor two. In addition, the 1-loop 
running is very different, since the 2-loop term in the /3-function contributes 25% and 
more in this range of r. We conclude that r ~ 0.17 fm is outside of the universal 
perturbative domain. 

We may also consider the perturbative relation between agg(r) and the finite volume 
coupling g^{L) that will be discussed further in the following section. This reads 



a Z + 0.9980 a^L) + a(L) 



9\L 



r)/(47r). Using the result of the 



following section, we get a, 



qq\ 



0.35i?o) = 0.39 + (9(a ), which is significantly below 
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the nonperturbative result agg(r = 0.35i?o) = 0.54(2). However, at this value of r, the 
0{a^) term is 30% of the 0(a) term and we arrive again at the conclusion that r is 
too large to use perturbation theory. Within our approach of extrapolating the lattice 
numbers to the continuum, we could not reach significantly smaller values of r. 

4. Setting the Physical Scale in the Computation of the Run- 
ning Coupling a(L) 

In ref. [6], a renormalised coupling g^{L), that runs with the box-size L was calculated 
for different length scales L} At the largest scale investigated in ref. [6], the coupling 
has the value g^{L) = 4.765. At this value one can make contact with the low energy 
scales of the theory in infinite volume: one calculates the product LE, with E some 
energy scale of the theory in infinite volume. The choice of ref. [6] was E = \fK ^ with 
K the string tension, since data for two small values of the lattice spacing existed for 
this quantity. This can be improved by using E = 1/R[c) instead. 

We describe now the computation of L/ R[c)^ extrapolated to the continuum limit. In 
a first step, we determine the bare coupling /3 as a function of the lattice size L/a for 
fixed g^{L) = 4.765. Then, we use R[c)/a at different /3-values in the same range, 
form the ratio L/R[c) and extrapolate it to the continuum a/L = 0. 



L/a 


/3 


l^impr 


5 


2.4971(23) 


2.5500(23) 


6 


2.5752(28) 


2.6082(30) 


7 


2.6376(20) 


2.6674(24) 


8 


2.6957(21) 


2.7155(21) 


10 


2.7824(22) 


2.7989(23) 


12 


2.8485(32) 


2.8697(39) 


14 


2.9102(62) 





Table 2: Bare coupling /3 from ref. [6] and bare coupling Pimpr for the 1-loop improved action vs. 
the lattice size at fixed g^(L) = 4.765. 

At the largest values of the coupling, a significant lattice spacing dependence of the 
step scaling function was found in ref. [6] (cf. fig. 2 in that reference). It was therefore 
to be expected, that L/R[c) shows lattice artifacts. We attempted to reduce these 
artifacts by using the 1-loop Symanzik improved action given in ref. [5]. We denote 
the results that are obtained using this action with a subscript "zmpr". The /3-values 

^For details on the definition and the computation of g^(L), we refer the reader to [5, 6]. 
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are listed in table 2 together with the ones without improvement that were already 
obtained in ref. [6]. 

Both /3(iv/a) |^2(2,)=4.765 and /3impr(-^/«) |p2(/^)=4,765 are almost linear functions of log(iv/a). 
Therefore, the (noninteger) values Lj a and Limpr/O' for the values of /3 where we know 
the force (table 3) are easily determined by interpolation. At these /3-values values, 
we calculated Rj;[c)/a as discussed in section 2.2. 

Examples of the data for Lj Rji^c) are plotted in fig. 4. It is well worth noting that 
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Figure 4: L/R^Jc) (o), Limpr / R^Jc) (•), L/R^Jc) (A) and Limpr / R^Jc) (filled A) as function of 
the lattice spacing a/L and c = 0.5, 1.0, 1.65 and 3.0 from top to bottom. The lines are the linear 
extrapolation to the continuum using a/Limpr < 1/5- At a/Limpr = the extrapolated values are 
plotted together with their error bars. 

fig. 4 is the first time that the success of the Symanzik improvement program can be 
demonstrated in 4-dimensional pure gauge theories: the dependence of the ratios on 
the lattice spacing is much weaker when the 1-loop improved action is used. In fact. 
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due to the weak cutoff dependence, the data with improvement can be extrapoiated 
to the continuum with confidence. 

The physicaiiy most interesting case is around c = 1.65 and we use i?(1.65) = i?o — 
0.5 fm in the foiiowing (at finite vaiues of the iattice spacing, Rj{c) depends on the 
orientation d. In the continuum limit, this dependence disappears and we drop the 
corresponding index.) Our extrapolation to the continuum yields 

L/Ro = 0.573(37) (4.1) 

or — using Rq = 0.49 fm as an illustration also for SU(2) Yang-Mills theory — L = 
0.281(18)fm. For comparison, we also give the numbers that were obtained from the 
force along the other two orientations: L/Rq = 0.587(39) for and L/Rq = 0.595(31) 
for J3. The latter extrapolations had to be performed without the /3 = 2.85 points. 

These extrapolations are stable: Within the error bars, we obtain the same results if we 
include one more point at larger a or if we remove the last point from the extrapolation. 
Furthermore, the difference of the extrapolated value from the point at largest a that 
was included in the extrapolation is only about 15%. Therefore, the above calculation 
gives the momentum-scale for the running coupling g^{L) to a precision of about 6% 
in the continuum limit. 

It is more easy to control systematic errors in the calculation of i?o than in the de- 
termination of the string tension. Nevertheless, it is interesting to check, in how far 
one obtains a different result if the scale is set by the string tension. We have there- 
fore determined the string tension from Fj;[r) = K -\- 7r/(12r^) for K > 1.8. This 
range in r was chosen such that the fi = 2.85 data could be included. The 7r/(12r^) 
correction[14] is less than 14% and the results of section 3.3 provide some evidence 
that the uncertainty on the correction may be neglected. Fxtrapolating Limpr^/K to 
the continuum as shown in fig. 5, gives 

L^/k = 0.675(60) (4.2) 

or (with \/K = 425MeV) L = 0.313(28)fm. The small change compared to the 
determination that uses L/Rq is due to the fact that the force in SU(2) Yang-Mills 
theory deviates somewhat from the Richardson model in the relevant r-range. 

The statistical uncertainty in both determinations of L in physical units is quite similar, 
because they originate from the force in about the same range of r. However, the string 
tension determination uses the assumption on the subleading correction to the force. 
Although fig. 2 provides reasonable evidence that we may use the correction derived 
from bosonic string vibrations, an uncertainty due to this assumption is not easy to 
quantify. 
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5. Conclusions 



It has already been pointed out in [6], that a good way to set the scale in the pure 
gauge theory is through the force at some finite distance well inside the nonperturbative 
region. We have shown here that this proposal can be carried out in practice and that 
the quantity Rq provides a low energy scale that can be calculated precisely. Compared 
to the string tension this avoids the extrapolation in r and is therefore independent of 
parametrisations of the force. The relation to experiment is based on the assumption 
of identifying the static QCD potential and the phenomenological potential. 

Using Rq and the 1-loop improved action for the simulation determining g^{L), we 
were able to compute the scale in the running coupling of [6, 5] to within 6% in the 
continuum limit. 

The force itself could be extrapolated to the continuum limit with all systematic errors 
taken into account. In order to do this, it was most important to have precise data at 
several rather small values of the lattice spacing. In addition, lattice artifacts could 
be suppressed by choosing a tree-level improved finite difference for the lattice force 
in terms of the lattice potential. 

i?o calculated in the way described here will clearly be a useful reference scale also in the 
SU(3) pure gauge theory. It should be mentioned, however, that once one approaches 
full QCD, there are other quantities that are directly measured in experiments that 
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can play this role. Nevertheless it will be interesting to reconstruct the force along the 
lines of section 3.2 and compare to phenomenological models. In addition, a check of 
the relations eq.'s (2.3,2.4) will be of interest. 
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him in particular for critically reading the manuscript. 
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A The Calculation of the Potential 

The force between static quarks is most efficiently calculated from "smeared" Wilson 
loops. We have applied a variation of smearing based on the original ideas of ref. [18, 
19]. The smearing and variational technique employed her deviates only slightly from 
the one used by Michael et al. [10, 7]. A detailed discussion of the merits of the 
method can be found in [7]. Here, we describe the exact implementation that was 
used and discuss the results. 

We start from lattice gauge fields U^{x) gSU(2), that are the gauge connections be- 
tween points X and x-\-jla^ // = 0, 1, 2, 3 of a hypercubic lattice A with lattice spacing a. 
In a first step, the space-components of the gauge field are "smeared" by one iteration 

3 

Uk{x) ^ U',{x) = V { Uk{x) + ^ [ U,{x)Uk{x + U)U]{x + ka) (Al) 

+ [/|(x — ja)Uk{x — ja)Uj{x — ja + ka) ] }, A; = 1, 2, 3, 

Here V denotes the projection into SU(2) and we have used io = 1/4. Next we block to 
a sub-lattice A^; characterized by a vector b with components bi = 0, 1. The sublattice 
has 1/2 the number points per space-dimension by restricting the components of x/a 
to be even [bi = 0) or odd: Xi/a = bi mod 2, z = 1,2,3 for x G A^;. The gauge fields 
on the sublattice are simply 

Uk{x) = Uk{x)Uk{x + ka), xeA^. (A2) 

On this blocked lattice we apply eq. (Al) Si times (with a ^ 2a and U — > U). In order 
to be able to calculate the potential along the three orientations di = 2a (1, 0, 0), = 
2a (1,1,0), = 2a (1,1,1), we then construct generalized gauge fields V^'^ that 
connect x = (xq, x) with ( ' + di), I = 1,2,3. yw are obtained by averaging the 
parallel transporters over the different shortest lattice paths on A^. 
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We computed smeared Wilson loops with space extent r = rii di, = 1, 2, Af using 
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Figure 6: \og[C{t,r)IC{t — a,r) at (] = 2.6. The data points are from top to bottom for r = 
5(i2, 4(i3, 6di and Srfs. The long dashes represents the fit to eq. (A6) and the full line is the fit to 
eq. (A7), plotted in their respective fit ranges. As an illustration, the fits to eq. (A7) are continued 
to t— values outside of the fit ranges as short dashed curve. 

the appropriate product of y^'-* for the space-like parts and the 1-link integral[20] for 
the time-like parts of the loops. Inserting V after smearing level Si at time and Sj 
at time t, one obtains a matrix correlation function. The latter corresponds to matrix 
elements [T^[f))ij^ n = t/a of powers of the transfer matrix T in the corresponding 
charged sector of the Hilbert space. ^ 

The generalized eigenvalue equation 

E[r'^+^(f)],.(f)f = A,(f) T.[Hr)hAr)f\ W) < W), ... (A3) 

i i 

^In this discussion we neglect corrections due to a finite time extent of the lattice and use the 
convention that the energy of the vacuum is zero. 
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gives estimates Xk{r) that converge, for n — )• oo, (exponentially) to the lowest eigen- 
values of the transfer matrix[21]. 

In our simulations (cf. section (3.1)), we found that the first two eigenvalues are fairly 
stable with respect to n, once n > 2. Higher eigenvalues are difficult to determine 
since - apparently - the different states generated by the above smearing procedure 
are not linearly independent to a sufficient degree. We used this variational technique 
to construct a correlation function 

C{t,r) = T.<rii'\HrlhAriP (A4) 

hi 

that obtains only small contributions from excited states because it is projected onto 
the approximate ground state. Furthermore, we estimated the gap in lattice units by 

a A(f) = log(Ai(f)/A2(f)). (A5) 



We then determined the potential from fits 

C{t,r) = C^{f)exp{-Vu{f)t) and (A6) 
C(t,r) = C^(f)exp(-V^(f)t)[l + ^exp(-A(r)t)] . (A7) 

Here A(r) was fixed to the value determined from the variational method eq. (A5). 
Vu{r) is an upper bound for the potential since the excited state contributions are 
neglected in eq. (A6). On the other hand, Vi{f) gives a lower bound for the potential 
provided A(r) as determined from eq. (A3) is a reasonable estimate for the gap. This 
is true, because eq. (A7) parametrizes all corrections to the lowest state by an effective 
contribution from the first excited state. It thus overestimates the corrections at large 
t (i.e. t-values outside of the fit range). As they are lower and upper bounds, one can 
- in principle - use the combination of eq. (A6,A7) for any fit ranges (in t), irrespective 
of whether the fits are statistically satisfactory. If this is done allowing too small values 
of t, however, Kt(r) and Vi{f) become significantly different and one has a dominant 
systematic error Kt(r) — Vi{f). In the present calculation, we have a relatively large 
range in t. We chose the fit ranges such that they yield a satisfactory (including the 
correlations of the data). This resulted in Kt(r) = Vi{f) within one standard deviation. 

We give examples of the data together with fits in the form of the time derivative of 
C(t,r) in fig. 6. The values of r chosen for the figure are in the range that is most 
relevant for the analysis of the force discussed in this work. It is evident from the figure 
that C(t,r) provides a correlation function which receives only small contributions 
from excited states in the accessible range of t. Therefore we can extract the ground 
state potential with moderate error bars. 
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Table 3: Force and improvement radius in lattice units. For convenience we have also included the 
results of [7]. In physical units r/ covers roughly the same range for all five fi. 
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